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1. INTBODDCTION AMD SIATBMEMT OP THE MAIM IBSULTS 

Th* "fiaultancout stablllMtlon problM** - la aitbar ditetata or 
conti&uoua tiae ~ conalata in anawaring the following quaatlont 


Given an r>tuple 6^(s)t...tG^(a) of p proper tranafar 
functions, does there exist a compensator K(a) such that the closed- 
loop systems Gj^(s)(I+K(s)Gj^(s)) ..,G^(8)(I + !C(s)G^(s)) ^ are 

(internally) stable? 


As pointed out in [13], this question arises in reliability theory, 
where G 2 (s),. ..,G^(s) represents a plant G^Cs) operating in various 
SK>des of failure and K(s) is a nonswitching stabilizing compensator. 

Of course, for the same reason, it is important in the stability analysis 
and design of a plant which can be switched into various operating modes. 
The simultaneous stabilization problem can also apply to the stabilization 
of a nonlinear system which has been linearized at several equilibria. 
Finally, it has been shorn [14], (20 1 that to solve the case r*2 is 
to solve the well-known problem considered by Youla et al in [21] : When 
can a single plant be stabilized by a stable compensator? This corres- 
pondence also serves to give some smasure of the relative depth of this 
problem. 

In order to describe the results obtained via this correspondence, 
we need some notation. First, set n^* McMillan degree of G^(s). In 
the scalar jj^pt^^output setting (m*p*l), we regard each G^(s) as a 
point in H ^ , viz. if 


GfCs) • Pj(s)/q^(s) , where 


n n -1 n 

Pj^(s) ■ a^j^ + . .. +a^ ^s , and qj^(s) •b^j + •••‘♦^b^^a +s 


“ni' 


th«n 6^(s) corresponds to ths vsetor 
Morsovsr, sines p. 








in ths open dsnss sst Rst(n^) c s 


and ^rs rslstivsly prins* this ssetor liss 
“ ^ (sss I 31 for ths strictly 

proper csss). In [14], Seeks and Murrey used ths techniques of frsctisnsl 
representations [8] snd the correspondence nentioned shove to give 
explicit inequalities defining the open set 


U c Rat(nj^) xRat(n 2 ) 

of pairs (G^Csl.G^Cs)) which are slaultaneously stabiliaable. In 
[20] Vldyasagar and Viswanadham shoiKd, using similar techniques, that 
provided max(n,p) > 1 the open set U of pairs (G 2 (s),G 2 (s)) which 
can be stabilized is in fact dense. 

This can be made precise by topologizing a point G^(s) in the set 
p-{pxm Gj^(s) ; degree Gj^(s)-nj^) 


(n +1) (iiq>) 

as a vector in B via its Hankel parameters: If 


Gj^(s) 


00 



then G^(s) corresponds to the n + 1 pxm block matrices n+1^ 

which determines G(s). It is known that ^ ^ is an (n(m*f p) ■fq>)-manifold 
(see [ 7 ], [12], [ 5 ]), although this is not important here. Hhat is 
Iflportant is that ^ is a topological space. 

*^,p 

One of our main results concerns the generic stabilizabllity problem; 
that is. 


Question 1.1 . Fix m,p,r, and n^. Is the set U of r-tuples 
G. (s),...,G (s) fdiich can be simultaneously stabilized open and dense 



« 


It is slso lai>ortsnt to ssk» for rsasons of global robustnsst of 
algorithms finding such a compensator « for eosg>ansators vlth a fixed 
degree of complexity. 


Question 1.2 . Fix m,p»r, and n^. What is the minimal value of q 
(if one exists) for uhich the set W^ of r-tuples which can be simul- 
taneously stabilized, by a coiq>en8ator of degree is open end dense 



It should be noted that, in the case r«l, Question 1.2 is an 
outstanding, unsolved, classical problem. In this paper, we prove: 


Theorem 1.1 . In either discrete or continuous time, a sufficient condi- 
tion for generic simultaneous stabilizability is 

max (m,p) ^ r (1.1) 

Indeed, if (1.1) holds, then the generic r-tuple can be stabilized by 
a compensator of degree less than or equal to q, %diere q satisfies: 

r 

q[max(m,p)+l - r] > ^ n -max(m,p) (1.2) 

i-1 

In the case r*l, it is unknown whether generic stabilizability 
implies generic pole-assignability; that is, whether or not these 
properties of m,n, and p are really different (see [4]). Perhaps 
not surprisingly then. Theorem 1.1 follows from: 
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Th«or«a 1»2 » A sufflclsnt condition for |«n«rlc •laultanoous polo-* 
aotlgnability it (1.1) » where the coopenaator K(a) can be taken to 
be of degree q eatitfying (1.2). 


Here, eiaultaneoue pole-aetlgnability neant the aaeignability of 
r sett of self-conjugate sets of nutters {e,.»...»t . .) ct. 

XI 

In fact sharper bounds on q can be obtained (see (18], (11]). Our 
proof relies on the recent pole-placeoent techniques derived for r«l 
by P.K. Stevens in his thesis (18], which contains an i^>roveasnt on 
existing results in the literature, see also [9], (17]. We shall prove 
Theorem 1.2 only in the strictly proper case; the proper case Involves 
more technical arguments from algebraic geometry which can be found in 
[11]. We shall, however, give an Independent proof of Theorem 1.1 in 
the nonstrictly proper case, based on the equivalence of generic stabili- 
sabillty and existence of a solution to a generic "deadbeat c<mtrol" 
problem, which we can solve if (1.1) is satisfied. This argument 
extends the argument given in [ 4 ] for the case r*l and q*0. 

Note that if r»l, then (1.1) is always satisfied in which case 
(1.2) implies 


Corollary 1.3 . (Brasch-Pearson [2]). The generic p^n plant G(s) 
of degree n can be stabilised by a compensator of order q, where q 
satisfies 


(q ’*-l)max(m,p) > n 


(1.3) 


If r-2 and max(m,p) > 1, then (1.1) is again satisfied, so we 
obtain rather easily: 


Corollary l.A . (VldyMagar-Viavanadhaa [201). If r-2 and Mx(n,p)>l» 
than tha ganaric pair (G^ (a) , 62 ( 0 )) la ainultanaoualy atabllliabla. 


Moraovar, in this caaa wa know an uppar bound on tha ordar of tha 
raquirad conpanaator. For axaapla, if ■<*p*2, r**2, then q can ha 
taken to aatiafy 


q k +n 2 • 2 

On the other hand, in [20] tha explicit conditions defining the cloaad 
aat 




-U 


of pairs not sisultanaously stabilisabla were derived. Such conditions 
can be derived from our proof, but instead we refer to [ 10 ]. where 
Theorem 1.1 (excepting (1.2)) is proved by interpolation methods also 
yielding a set of explicit conditions in the range r < amx(m.p). 

Finally, we prove that the condition (1.1) is sharp in the following 
sense. 


Theorem 1.5 . If min(m.p)«l. then for fixed m.p.r and n^ the following 
statements are equivalent for proper plants: 

r 

(i) q€W satisfies q(max(m.p ) + 1 - r) + n . ; 

i -1 

(ii) the generic r- tuple 0 ^( 0 ) .. . . .G^) is simultanemtsly 

stabilisable in discrete or continuous time by a compensator 
of degree f q: 

(iii) the generic r-tuple Gj^(s) G^(s) is ainultansoualy 

stabilisable in discrete or continuous time. 
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In the strictly proper case It follows that (l)**(ili) is also 
equivalent to generic siaultsneous pole assignability. This bolds in 
the proper case as well, but requires a separate arguaent [113 • 


Corollary 1.6 . If nln(n,p)«l and r < bsx(b,p) then the generic 
r-tuple is simultaneously stabilizable by a coopensator of order 
precisely given by the least integer q satisfying (1.2). 


As a further corollary, we obtain one of the results obtained by 
Seeks and Murray in [ ], see also [15]: 


Corollary 1.7 . (Saeks-Murray) . Suppose m"p*l and r«2. Simultaneous 
stabllizability is not a generic property. 


We remark that these results hold also over the field € of 
complex nu9d>er8 - in particular, the couq>lex analogue of Corollary 1.7 
dispels a folklore conjecture concerning simultaneous stabilization 
using compensators with complex coefficients. 

Finally, over any field, the method of proof of Theorem 1.2 gives 
linear equations for a compensator simultaneously placing r(n<fq) poles 
when the generic hypothesis is satisfied. 


2. POLE PLACEMENT AMD THE GENERALIZED SYLVESTOR MATRIX: A PROOF OF 

THEOREM 1.2 

In this section we proceed to prove Theorem 1.2. Note that 
Theorem 1.1 and Corollaries 1.3 and 1.4 follow Imnediately in the 
strictly proper case from this theorem. Without any loss of generality 
we can assume that m b p, for, if R(s) stabilizes G^(s) than K^(s) 


omelNM. PAGE ® 


stabilizes G^(s). 

Suppose, first of ell, that p**!, so that ve are given a set of 
r, B input 1 output plants of McMillan degree ( n represented as 
f k i ? k i r k 1 

i«0 , i«0 , , 1"0 




i-O 


'»fpl 




( 2 . 1 ) 


for k ■-1,2, . . . , A. A 1 input, m output, compensator of McMillan degree 
S q is represented as 


fa.,.* f a,,.* f , ,.* 

1=0^* , lip tip »*^~*-* 


lip 




iip-i>* 


i-0 ®^i® 


( 2 . 2 ) 


Note that in (2.1) and (2.2) the coefficients p^^ V k and a^^ 
has been defined up to a nonzero scale factor. Moreover, for a strictly 
proper plant or compensator, p^ ■ 0, a. • 0 V J ■ 1, . . . ,m-»-p - 1; k ■ 1, . 

The associated return difference equation, det(I -<-K(s)G^(s) ■ 0 
is given by 




(2.3) 


V k - 1,2 r 


A generic r-tup.l.e of plants define a mapping X> via equation (2.3), 
between the plant parameters and the coefficient of the return difference 
polynomials given by 


^ . jj(q+l) (m+p) ^ j^r(n+q+l) 


(2. A) 


where 
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P 


0 


X (Aq A ) - (A0.....A ) 


PO 


(2.5) 


where 


A, - , a2^ , , a^^) 


( 2 . 6 ) 


"l- 


1 

1 

r "" 

Pli 

Pli* • 

• • Pli 

^1 

2 

r 

P2i 

• • • e 

P2i* • 
• • • • 

• • P2i 

• • • • 
_1 

e • • • 

2 

r 

^m+pi 

Pm+pi’ 

P|B+p^ 


(2.7) 


The matrix in the right hand aide of (2.5) is classically known as 
the generalized Sylvester matrix and Is of order (q-fl) (m4-p) x r(n-fq*fl). 
For r«l its rank has been analyzed by Bitmead, Kallath, Kung In [ 1 ]. 

In particular, for a generic plant, it is known to have full rank. For 
rkl, we have the following: 


Lemma 2.1 . The generalized Sylvester matrix is of full rank for a generic 
r-tuple. 


Proof : See Appendix I. 
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Tamm AttuM ■ln(tt»p) • 1» A suffieiMt condition for gmorie 

polo «oolgnmont» for on r tuplo of otrietly propor plonto by o propor 
conptnsotor io givon by 

r 

(q + l)(n + p-r) b Jn.-r + l (2.8) 

i-1 ^ 


Proof ; Ue prove this Lonno aoouaing for notationol convoriancc that 
n^"n V l*l,...,r and analyse the napping x n* defined by (2.A}« 
(2.5). Assune 




k 

**n+p,n 


1 V k-1 


r 


and that the coefficient of in all the r return difference 

polynomials (2.3) has been normalised to 1. 

Thus a sufficient condition for generic pole aseignnent is that X* 
is onto. Here the mapping 

( 2 . 

is given by 


X (Aq, . . . ,A^_^A* ) ■ (Ap *q-l^q^ 


PO 


e • • a 


n I 

I 

*U 


(2.10) 


where 


K “ ^•lq‘2q’---»*m^-lq> 


and p* is obtained from p^ by deleting its (at+p)*^ 


row. 


n 
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By Lcosm 2.1 the Mtrlx in the right-hand aida of (2.10) ia of 
full rank for a ganaric r-tupla of planta* and has tha ordar 
(q^DCa-fp-l) X r(n4q). Tharafort* a aufficiant eondltioa for 
ganaric pole placanent ia given by 

(q + l)(a4-p) - 1 k r(n + q) (2.11) 

which is sane as (2.8) for n^^^n V i*l«...»r. 


The proof of Theoren 1.2 now proceeds by a reduction to the case 
Bin(a,p)-1, which has been traated in Laanas 2. 1-2. 2. This procadurs, 
which is called "vectoring down", is adopted from the case r-1, 
studied in P.K. Stevens thesis [IS]. 


LeiMMi 2.3 . Given an r-tuple of pxn plants G^(s) of degrees n^^, 
each with distinct simple poles, there is an open danse set of ixp 
vectors v € such that vG^(s) has degree n^. 


Proof; If r»l, then we nay expand G(s) 


n K. 

=<•> • I rrp 

i-1 • 1 

in a partial fraction expansion, where € t and each has rank 1. 
Now, the set of real vectors v such that vR^^ ^0 is clearly open 

and dense in . Defining U.,...,U similarly, set 

A n 

n 

V* r\^ 4 ‘ Thus, V is an open dense set of vectors with the required 
1-1 ^ 
property. 

If r>l, one obtains, ss above, sets Vj^,...,V^ in Wp having 


M open dense intersection ^ V. . 

i-1 ^ 


Q eC a Ue 
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!•■■■ 2.4 . Civtn an r-tupla of pxa plants G^(a) thara axista a 
constant gain output feedback k such that the closed loop systsM 
C^(s)(1 4kG^(s))~^ have distinct single poles. 


Proof I For r«l. the set of K such that the closed loop systea 

has sinple poles is the cot^leaent in of an algebraic set. It 

is well known [ 2 ] that this set is nona^ty; therefore, W. is open 

r 

and dense. Taking any K in the open dense set gives the 

1-1 ^ 

desired conclusion. Q.E.D 


Thus, choosing any (v,K) C 1R^ x we have a napping fron an 
an open dense set 

l,^ * I" X ... X ■*’ I" 1 ^ I!! 1 

(v,k) ^m,p ^n,p ^«,1 ^,1 


♦(.,»)“=! ''’’ I .! • (vCj<.)(«.KGj(.»-^J.i 


which is rational in the Hankel parasietcrs 

Leases 2.1-2. 2 to the case siin(st,p) - 1, i.c. 10 1*^ 

gives - via conposition with ^ - an open dense set of 

r >'• 

^n.P 


(G^> . 


Applying 


y" which can be sinultaneously pole-assigned. 
‘T»,P 


Q.E.D 


3. GENERIC STABILIZABILITY CONDITION OF AN r-TUPLE OF PROPER PLANTS 

In this section we proceed to prove Theoren 1.1 independent of 
Theoren 1.2. Ve first show that the following three stateaents arc 
equivalent . 
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I. A generic r-tuple of proper plante le eteblllsable with reepect » 

to the open left helf plane by a proper eoapaneator of degrM ( 

II. A generic r*tuple of proper planta la atablliaable with raapect to 
the interior of the unit diac» by a proper coapeneator of degree d q. 

III. A generic r**tuple of proper plante la pole aasignable at the 
origin by a proper coapeneator of degree d q. 


Leaatf 3.1. I <■> 11 


Proof ; Conaider the conforaal tranafomacion 

♦(a)- (a + l)/(a-l) (3.1) 

which aapa the r-tuple of proper planta gj^»g2 **«*»tf the r-tuple 

of proper planta g^*****!^ where g|(a) ■ g^(^(a)) except for the 
algebraic aet of plante aatiafying - *'g^(a) hea a pole at a - 1 for 
aoae i-l,...tr". The proof now followa froa the two facta. 

1. ♦(a) aapa the open left half plane onto the interior of the unit 
diac. 

2. The aapping 

(Sj*****!^) 

and ita inverae, aa; the generic r-tuple of proper planta to the 
generic r-tuplc of proper plante. 
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t^imna *1.2. II <•> 111 


Proof ; Sufficiency is clear and follows by an analagous argument of 
Lemma 3.1 with <|)(8) »s + a, a>0, a€lR. 

To prove necessity, we have the following: For each r*l,2,... 

(shown easily by assuming statement II and considering a > 0, 

a € IR) . 3 an open dense set of of r-tuple of plants for which 

there exist a compensator of degree <. q which places the poles in the 
interior of the disc of radius 1/r centered at the origin. 

Consider the set 


OO 

U - 0 

r-1 

Clearly, U is a dense set by the Baire Category Theorem [13]. Since 
the mapping x given by (2.4 is linear, it has a closed image. Moreover, 
every r-tuple of plants in U admits a sequence of compensators which 
places the poles arbitrary close to the origin. Since image of X Is 
closed, U is contained in a set V of all r-tuple of plants for which 
there exists a compensator which places the poles at the origin. By 
the Tarski [19] - Seidenberg [16] theory of elimination over IR, V is 
indeed defined by union and/or intersection of sets given by polynomial 
equations or Inequations f^ > 0, fg“0. Finally, since U is dense 
in V, fg(U) * 0 ■> f g - 0 so that V is defined by strict polynomial 
Inequalities. Hence V is open. Moreover, since U is dense, V is 
also dense. 


l^nmia 3.3 . For a generic r-tuple (r^m + p) of min(m,p) * 1 plants 


III <*> (q + l(m + p-r) > r(n-l)+l 


Proof ; The only nontrivial part is to prove sufficiency for the case 

r(n + q) < (q + l)(m+p) < r(n + q + l) 

(The other cases follow easily from the fact that the associated 
Sylvester's matrix is of full rank for a generic r-tuple.) 

To prove sufficiency, for the above case we want to show that the 
vector 


( 0 , 0, *••*•••*, 0 , 8 ^, 82 , 

■4— r(n +q) ■ r ■» 

indeed belongs to the image of X (defined by (2.5)) for some 
8^ 0, 1 ■ 1, . . . ,r. 

Partition the Sylvester's matrix in (2.5 as [Sj^S 2 l where S^^ is 
of order (q + l)(m + p) x r(n+q). Clearly we are solving the pair of 
equations 


[A_, . . . ,A ]S. * [0,...,0] (3»2) 

u q 

t Aq , . . . , ]S2*[Sj^,...,s^] (3.3) 

We claim that for a generic r-tuple of plants (3.2) has a solution for 
a nonzero vector A for otherwise if A * 0 we have 

q q 

[Aq A^_j^]S’ - (0,...,0) (3. A) 

where S| is of order q(m + P) ><r(n + q) obtained by deleting the last 
m + p rows of Sj^. From (3.4) (A^ ^q-1^ * — fg of full 

rank genetically. Thus the only solution of (3 2 ) is the zero vector 
which is a contradiction since the kernel of S^ is at least of dimension 
1. On the other hand, for A^ j* 0, for a generic r-tuple of plants 
the right-hand side of (3.3) is a vector none of whose entries are zero. 



lb - 


Theorem 1>1 then follows from LeanB 3.1» 3.2» 3.3 end the 
vectoring down technique used in the proof of Theorem 1.2 in 
Section 2. 


4. PROOF OF THEOREM 1.5 

To say there exists q € K satisfying (1.2) is to aiax(m,p) > 
Thus, (ii) follows from (i) by Theorem 1.1. 


(ii) (ill) since (111> Is weaker than (ii). 

By Lemma 3.1, in order to prove (ill) •> (1) it suffices to assume 
that G^(s) ,. . . ,G^(s) are simultaneously stabilizable in continuous-* 
time. 


Proposition 4.1 . The generic (ro + D-tuple of 1 xm proper continuous 
time plants of degree n is not simultaneously stabilizable by a proper 
compensator of finite (but not a priori bounded) degree. 


Proof: Consider the domain of (simultaneous) stability 


n^+q 


•a. • 

^ e,c ) : y C. .s^ has all roots in D, } 

11 In xo i 




n -fq 
r ^ 


and Its convex hull c IR x ... x IR . Clearly » a necessary 

condition for generic simultaneous stabllizabllity is 


lmage( X^) H ^ 0, 


for an open dense set of n* Since 
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n(2>) C {(C^^) ; >0} 


it will suffice to prove: 


Lemma 4.2 . If r“m + p, then there exists an open set of r-tuples 
n such that image ( X^) contains no vector with only positive entries* 

We fix the value of q and construct the associated Sylvester 
matrix S. We claim that the open set E of plants defined by 


E - {(P«,P ,P_)|Pn‘ V j*l,...,n has all the entries negative} 

0 1 n 0 3 

cannot be stabilized by a proper conq>ensator of degree i q. 

Suppose the above is not true, then there exist n€E, such that 


image ( x^) ^ 

or in other words V i*l, r(n + q + l) and 

a S ■ a 


has a solution. Writing S as 

S - [S' S"1 

where 



and Pj “ 0 for all j > n 


(4.1) 


(4.2) 
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Equation (A.l) can be written as 

a’ II I S’~^ S”1 - a (4.3) 

where S' ^ is given as follows 



The identity matrix of order (q + l)(m + p) in (4.3) forces a' to have 
all the entries positive. Moreover, since n€E,S* ^ S" has all its 
entries negative so that a'(S'“^ S") has all the entries negative which 
is a contradiction since a is a positive vector. 

Finally it is shown that E is not an empty set. For a fixed 
Pq • P* choose the vector 6 to be so that Pq ^6 has all its entries 
negative . Let 

P* • (6,6.. ...6) j - l,...,n 

■*— m + p -■ > 



OF POOR 

80 that 

(Pq.Pi ^ ® 

Q.E.D. 

Ranark ; If image ( x.^) affine hyperplane, then the necessary condition 

imageC X^) i 0 

of course is sufficient, i.e. implies 

image( Xq) ¥ 0 

This fact was used by Chen, together with 

Lemma 4.3. (Chen [6]) If r“l, fiO*)“{(c, ,...,c_) : c^>0} 

1 n 1 


PRG^ 'S 

qUAUT( 


to give precise conditions for stabllizability in the case 
r*l, q ■ 0, min(m,p)*l, and tnax(m,p) • n - 1. This technique can be 
adapted in the cases r > 1 to give explicit conditions - in certain 
cases - defining the open set of simultaneously stabilizable plants 
when r>max(m,p), see [11]. 

Note that Corollary 1.6 now follows from our previous results on 
the generic rank of the generalized Sylvester matrix, while Corollary 1.7 
follows either from Theorem 1.5 or Proposition 4.1. 
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APPENDIX I: PROOF OF LEMMA 2.1 

TIte generalized Sylveator aatrix is co-ordinatized by r(n-f l)(n-f p) 
parameters, and it is sufficient to show the existence of one principal 
minor with nonvanishing determinant. 

By reordering the rows and columns, the generalized Sylvestor matrix 
can be written as 


where 


S - tQi»Q2*‘ “’V+p^^ 
^i * ^^11*^12’*’* *^ir^ 


( 1 ) 

( 2 ) 


jk 


JO ‘^jl 


■Jn 


k • k 
’jO *^J1 


Jn 


(q + 1) 


(n + q + D- 


(3) 


in the notation of (2.7). Moreover, each p^j^ is referred to as a 
•block’ of S. 


Define a set M of matrices as follows: "n belongs to M 

provided m is obtainable from one of the matrices p^^ in (3) 
either by deleting the first columns or the last Q 2 rows 

> 0 ." 


Proposition A.l . 


exists a principal minor 

* 

* *^jn 


Every element m of M 

m 


m e M of 
P 


has the property that there 
* 

, a coordinate p and an Integer 


j such that p 
''m 


is a summand in det m where j 

p tn 


is the order 


of the minor. 
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Proof ; 


Clear from the structure of 



The following Is as algorithm to construct a principal minor 
with nonidentically-vanlshlng determinant. 


Algorithm: 


Set S ■ S , Initialize ^ * 0 


1. Set 5-C+l- 


2. Look at Obtain the principal minor m of P^^^. satisfying 

Proposition 1. If there is more than one possible choice, choose 

A it 

the one containing the first column. Define ctr“p and Jr"j • 

ID S ID 

* 

3« Delete the rows and columns corresponding to the coordinate p from 

m 

S. Renumber the blocks of the resulting matrix and call it S. 

(Every block of S is to be identified as a minor of the correspond- 
ing block in S obtained by row or column deletion.) 

4. Do the sane "delete" operation as in step 3 in S. 


5. If S is empty, terminate. Otherwise go to 6. 


6. Set k«C- 

Construct the principal minor m^ of S by choosing those elements 
of S whose corresponding row and column has been deleted in Step 6. 


Proposition A, 2 . During the execution of the above algorithm, S can 
always be decomposed Into blocks belonging to M. 


Proof: Clearly S satisfies the above proposition, since each block 

p belongs to M. Each Iteration of the algorithm deletes either the 

J ^ 
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first columns of the first block column of S or the lest 02 
rows of the first block row of S. The proposition thus follows from 
the definition of M. 


Proposition A. 3 . constructed in Step 6 of the algorithm has a 

nonldentically-vanishing determinant . 


Proof ; Wfe prove the proposition by showing that det m has a sumsMnd 

k ^ 

given by Ha, in the notation of the algorithm. This is clear, 

1-1 

however, by observing that in the C iteration the matrix S has a 
principal minor, the determinant of which has the summand II a. , 

i-5 ^ 


where k is defined in Step 6 of the algorithm. 


, - T 


4 
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